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CO 
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The usual formulas for the correlation functions in orthogonal and symplectic ma- 
trix models express them as quaternion determinants. Prom this representation one 
can deduce formulas for spacing probabilities in terms of Fredholm determinants of 
matrix-valued kernels. The derivations of the various formulas are somewhat involved. 
In this article we present a direct approach which leads immediately to scalar kernels 
for the unitary ensembles and matrix kernels for the orthogonal and symplectic ensem- 
bles, and the representations of the correlation functions, cluster functions and spacing 
distributions in terms of them. 



1. Introduction 

In the most common models of random matrices the eigenvalue distribution is given by 
a probability density P^{x\, ■ ■ ■ ,xn). If F is a symmetric function of iV variables and the 
eigenvalues are Ai, • • • , Ajv then the expected value of F(Xi, ■ ■ ■ , Ajv) is given by the formula 



J ■■■ J F(xt, ■ • ■ ,x N ) P N {x u ■ ■ ■ ,x N ) dxi 



■ dx 



N- 



Here dx denotes Lebesgue measure on a set which necessarily contains the eigenvalues (R 
if the matrices are Hermitian, T if the matrices are unitary, C if the matrices have gen- 
eral complex entries). The function P N (xi, ■ • • ,x^) gives the probability density that the 
eigenvalues lie in infinitesimal neighborhoods of x\, • • • , xpj. The n-point correlation function 
R n (xi, - ■ ■ , x n ) is defined by 

Rn{Xi, • • ■ , X n ) = — J " ' J ^N\ x li ' ' ' ? x n) dx n+ i ■ ■ ■ dx??. 
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It is, loosely speaking, the probability density that n of the eigenvalues, irrespective of order, 
lie in infinitesimal neighborhoods of xi, ■ ■ ■ , x n . (It is not a probability density in the strict 
sense since its total integral equals N\/(N — n)\ rather than 1.) 

Sometimes easier to compute are the n-point cluster functions T n (xi,---, x n ). These are 
defined in terms of the R m (xi, • • • , x m ) with m < n. Reciprocally, the R n (xi, - • ■ , x n ) may 
be recovered from the T m (xi, ■ • • , x m ) with m < n. These will be discussed in Section 2. 

Fundamental to the study of spacings between consecutive eigenvalues is the quantity 
E(0; J) (the E here does not represent expected value), which is the probability that the set 
J contain no eigenvalues. More generally E(ni, ■ • ■ , n m ; J 1; • ■ • , J m ) denotes the probability 
that each J, contains precisely n« eigenvalues. A related quantity is Pj(xi, ■ ■ • , x n ), the 
probability density that the eigenvalues contained in J lie in infinitesimal intervals about 
xi, ■ • • , x n . These are useful for the study of spacings between several eigenvalues. 

All the functions R n and T n , as well as the just-mentioned probabilities, are obviously 
expressible in terms of the density function P^. But since we generally think of iV as 
large (indeed we are often interested in "scaling limits" as iV — > oo) it is desirable to find 
alternative expressions for these quantities which do not increase in complexity with large 
N, expressions in which N appears simply as a parameter. 

These have been found for some particularly important matrix ensembles. They are 
ensembles of Hermitian matrices where P/v has the form 

P N (Xx, - ■ ■ ,X N ) =C N Y[ \Xj~X k f Y[w(xj), 
j<k j 

where (3 equals 1, 2, or 4, w(x) is a weight function and cn is the normalization constant 
required to make the total integral of P/v equal to one; and ensembles of unitary matrices 
where P^ has the form 

j<k 

and the integrals are taken over any interval of length 2ir. The latter are Dyson's circular 
ensembles. (Of course one could also introduce a weight function here.) The terminology 
here can be confusing. The (3 = 1,2 and 4 ensembles are called orthogonal, unitary and 
symplectic ensembles, respectively, because the underlying measures are invariant under 
actions of these groups. In the Hermitian case the (3 = 1 ensembles consist of real symmetric 
matrices and the (3 = 4 ensembles consist of self-dual Hermitian matrices. For a discussion 
of these matters and the analogues for the circular ensembles see ||. 

Here is a rough outline of how the desired expressions were obtained when (3 = 2, which 
are the simplest cases. For Hermitian ensembles 

P N {xi, ■■■,x N ) = C N Y[(Xj - x k ) 2 Y[w(xj). 

3<k j 

On the right side we see the square of a Vandermonde determinant (for the circular ensemble 
it is a product of Vandermonde determinants), which is equal to the determinant of the 
Vandermonde matrix times its transpose. Performing row and column operations and writing 
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out the product, one obtains a representation 



P N (xi, x N ) = — det {K N (xj, x k )) j)k=1) ... tN (1.1) 

with certain kernels Kn{x,h) which are expressed in terms of the polynomials ortho normal 
with respect to the weight function w(x). These kernels are shown to satisfy a basic family 
of integral identities which are used to deduce that 

Rn(xi, ■ ■ ■ , x n ) = det {K N (xj, x k ))j,k=i,-, n - (I- 2 ) 

By manipulating the iV-fold integral that gives E(0; J) and using these expressions for the 
R n one deduces that E(0; J) equals the Fredholm determinant of the kernel Kn(x, y) acting 
on J. One obtains expressions for the more general quantities E(ni, ■ • • , n m ; Ji, • • • , J m ) 
in terms of the same kernel. Observe that iV appears only as a parameter in the kernel 
K N (x,y). 

For (3 equal to 1 or 4 the situation is less simple. From the point of view of the above 
outline, the problem is to concoct a kernel satisfying (|1.1|) and the family of integral identities 
which led to formula (|1.2j ). This problem was solved by Dyson || in the case of the circular 
ensembles and required the introduction of a new concept — the quaternion determinant. A 
quaternion- valued kernel was produced for which (|1.1|) and (1.2) held if all determinants 



were interpreted as quaternion determinants. Subsequently it was shown by Mehta 0] how 
to obtain analogous quaternion determinant representations for the Gaussian ensembles, the 
ensembles of Hermitian matrices with weight function w(x) = e~ x , using skew-orthogonal 
polynomials. This was generalized to general weights by Mehta and Mahoux ||. (See also 
||, [|, 00 After these quaternion determinant representations for the R n it is possible 
(by an argument to be found in A. 7 of 0, although not explicitly stated there) to deduce 
representations for E(0;J) 2 as Fredholm determinants of 2 x 2 matrix kernels These 
matrix kernels are the matrix representations of the quaternion kernels. 

If all one is interested in is the spacing distributions then this is certainly a very round- 
about way of obtaining Fredholm determinant representations. The purpose of this article 
is to present a direct approach, one which leads immediately to the scalar kernels when 
(3 = 2 and the matrix kernels when (3=1 and 4, and the representations of the correla- 
tion functions, cluster functions and spacing distributions in terms of them. It uses neither 
quaternion determinants nor a family of integral identities for the kernels. (That the correla- 
tion functions are equal to quaternion determinants for the (3=1 and 4 ensembles becomes 
a consequence of the representations.) 

What we do use are the following three identities which represent certain iV-fold integrals 
with determinant entries in terms of N x iV or 2N x 2N determinants with integral entries: 

det(0 i (x fc )) iifc= i i ... i AT ■ det(-0j(^it))j,ife=i,...,Arc/a:i ■ ■ -dx N 

= N\ det ( I <Pj(x) ip k {x) dx) . k=i n . (1.3) 
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det((j)j(x k )) jtk=li ... jN dxi--- dx N 
Pf(/ J sgn(y - x)fa(x)<l> k (y)dydx) . (1.4) 



det((pj(x k ) ipj{x k )) j=1; ... !2N , k =i,-,N dxf- dx 



N 



(2N)\ Pf ( J (hix) m%) - M%) M*)) dx ) j , k=1 ,., m - ( L5 ) 



Here "Pf" denotes Pfaffian (its square is the determinant) and ( |1.4|) holds for even N and 
must be modified for odd N. These hold for general measure spaces; in ( |1.4j ) the space must 
be ordered. Identities (|1.4j) and (|1.5| ) are due to de Bruijn ||, who traces (|T73|) as far back 
as 1883 0. 

In the application of ( |1.3| ) to the (3 = 2 case we are led to a matrix whose j, k entry equals 
5j ik plus the integral of a product of a function of j and a function of k. In our application 
of ( |1.4| ) and ( |1.5| ) to (3 — 1 and 4 the integrand is a sum of two such products. This, in 
a nutshell, is why 2x2 matrix kernels arise. What we do in each case is use one of the 
identities to express 

j ■■■ j Pn{x\i ■■■,xn) n( 1 + f( x i)) dx ^ ' --d^N (1-6) 

as a determinant or Pfaffian whose entries are given by one-dimensional integrals. Manipu- 
lating the integrals leads to a scalar or matrix kernel K^(x,y) such that the above iV-fold 
integral ((3 = 2) or its square (/3 = 1 or 4) is equal to det (/ + K N f). Here K N denotes 
the operator with kernel K^{x,y) and / denotes multiplication by that function. Taking 
/ = — Xj gives immediately the representation of E(0; J) or its square as det (I — K^ Xj). We 
also obtain representations for the more general quantities E(ni, • ■ ■ , n m ; Ji, ■ • • , J m ). Taking 
/ to be a linear combination of delta functions leads to representations of the correlation 
and cluster functions in terms of the matrices (K(xi,Xj)). These are 2x2 block matrices if 
(3 = 1 or 4. 

All the derivations follow the same pattern. First we write the integrand in ( |1.6|) in terms 
of a determinant or product of determinants. These representation are exactly as can be 
found in or 0, for example. Then we apply whichever of ( |1.3| )-( |i~5D is appropriate, and 
obtain a determinant whose entries are integrals. Finally, after writing the integrand as a 
matrix product if necessary, we use a general identity to express the last determinant as an 
operator determinant. 

The article is organized as follows. In section 2 we define the cluster functions and dis- 
cuss their relationship with the correlation functions. In sections 3 and 4 we show how the 
operator determinant representations of the integrals ( |1.6|) lead to formulas for the correla- 
tion and cluster functions as well as the spacing probabilities. In the following sections we 
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derive these determinant representations for the (3 = 2 ensembles, the (3 = 4 and 1 circular 
ensembles, and then the (3 = 4 and 1 Hermitian ensembles. We do things in this order since 
the circular ensembles are simpler than the Hermitian (because the exponentials are eigen- 
functions of the differentiation operator) and (3 = 4 is simpler than (3 = 1 (because there is 
a single integral in ( |1.5| ) and a double integral in ( |l-4j )). In the final section we derive the 
matrix kernels for the (3 = 4 and 1 Gaussian ensembles, which was the starting point of 

We emphasize that there are almost no new results in this article — perhaps the determi- 
nant formula ( |4.2| ) for the spacing probabilities and the more general form for the matrix 
kernels obtained in sees. 8 and 9. It is the methods used to derive them which are new and, 
we believe, show how they form a coherent whole. 

2. Correlation and Cluster Functions 

The n-point cluster functions T n (xi, • • • , x n ) are defined as follows. For each non-empty 
subset S of {1, • • • , iV} write 

R s = Rn{xiv- ■ ■ ,Xi n ) if S = {h, ■ ■ -,i n }. 

Then 

T n (x lt •••,*„)=£ (-l) n " m (m - 1)1 R Sl ■ ■ ■ R Sm , (2.1) 

where the sum runs over all m > 1 and, for each m, over all partitions of {1, • • - ,n} into 
nonempty subsets S\, - ■ ■ , S m . If we know the T n then the R n may be recovered by use of 
the reciprocal formula 

R n {x 1 , ■ ■ ■ , x n ) = £ {-l) n - m T Sl ... T Sm . (2.2) 

This is most easily seen by looking at the relation between generating functions for these 
quantities. 

Let A be a formal power series in variables Z\, • • • , Zn without constant term. For each 
nonempty subset S of {1, • • • , N} let A s denote the coefficient of Tiies z i m A. Let 

oo 

f( Z )= 

m=l 

be a formal power series in a single variable z without constant term, and set B = f(A). 
For each S we have 

Bs = y Zm\f m A Sl ---A Sn , (2.3) 

where the sum is taken over all m > 1 and, for each m, all partitions of S into nonempty 
subsets 5*i, ■ • • , S m . If fi ^ then f(z) has a formal inverse g(z). Since A = g(B) it follows 
from the above that for each S 

A s = Y J m\g m B Sl .-.B Sm . (2.4) 

We apply these fact to f(z) = — log(l + z), g(z) = e~ z — 1, taking any formal power 
series A such that As = R n {xi i: ■ ■ ■ , Xi n ) when S = • • • , i n }. Since f m = (— l) m /m in this 
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case, we see by definition ( |2.1| ) that with the same S we have Bs = (— l) n T n (xi 1? • • • , Xi n ). 
Since g m = (— l) m /m! (|2.4j) gives the reciprocal formula (|2.2| ). 

Observe that R n (yi, ■ ■ ■ , y n ) equals the coefficient of Z\---z n in the expansion about 
Z\ — ■ ■ • — z n — of 

/ ■ ■ ■ / Pn(xi, - ■ ■ , xtv) ]J[1 + ^ z r - y r )] dxi • • ■ dx N . (2.5) 

i r=l 

Here (1 + £)" =1 z r $( x ~ Vr)) dx denotes Lebesgue measure plus masses z r at the points y r . 
Thus we may take A(zi, • ■ ■ , z n ) to be this integral minus 1. Since B = — log(l + A) we see, 
after recalling the relationship Ts = (— l)^Bg, that T n (yi, ■ ■ ■ ,y n ) equals (— times the 
coefficient of z\ ■ ■ ■ z n in the expansion about z\ — ■ • ■ — z n — of 

log(/---/ P N (x u ---,x N ) Y[[i- + J2z r S(xj-y r )] dxx-'-dx^. (2.6) 

j r=l 



3. Formulas for the correlation and cluster functions 



Suppose we have proved that for each (3 = 2 ensemble there is a kernel K]y(x,y) such 
that for general / 




Pn(xi, • • • , x N ) JJ(1 + f(xj)) dx\ ■ ■ ■ dx 7v = det (/ + K N f). (3.1) 



As before Kn denotes the operator on L 2 with kernel Kn(x, y) and / as an operator denotes 
multiplication by that function. 

Take f(x) = J2r=i z r $( x ~ Vr) first. Then the operator / + K N f becomes the matrix 

[8 rs + K N (y r , y s ) z s ) . (3.2) 

V / r,s=l,—,n 

This may be seen by passing to a limit or observing (as we shall) that this is what happens 
during the derivation. It is easy to see that the coefficient of z\ ■ • • z n in the expansion of the 
determinant of fl3T2] ) equals det (^K N (y r ,y s )^j, and so by fl2.5p we obtain formula ( |1.2| ) for the 
correlation functions. 

For the cluster functions we use (|2.6|) and the general expansion 

oo t i \m+l 

log det (I + K) = tr log (I + K)=Y tr K m 

m 

to deduce 

T n (yi, ■ ■ ■ , y n ) = - J2 K N (y a i,y a2 ) ■ ■ ■ K N (y an , y al ), 
the sum taken over all permutations a of {1, • • • , n}. 
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Turning now to the (3 = 1 and 4 ensembles, suppose we can show that for each of them 
there is a matrix kernel Kjsi{x,y) such that 



J ■■■ J Pn{xi,---,xn) n( 1 + f(x j ))dx 1 ---dx N = ^det (I + K N f) . (3.3) 

3 

Then taking f(x) = J2r=i z r ^{ x ~Vr) we see by (|2.5|) that R n (yi, • • • , y n ) equals the coefficient 
of z\ • • • z n in the expansion of 



det [5 rjS + K N (y r ,y s ) z s ) . (3.4) 

V / r,s=l,---,n 



This coefficient is not nearly as simple as in the (3 = 2 case. However (|2.6|) shows that 
T n {yi, • • • , y n ) equals (— l) n+1 times the coefficient of z\ ■ ■ • z n in the expansion of 

i log det (5 r>s + K N (y r , y s ) z r ) . 

Thus we obtain 

T n (yi, ■ ■ ■ ,y n ) = — Y, tT K N{y*i,ya2) ■ ■ -K N (y an ,y al ), 

which is hardly more complicated than for (3 = 2. To obtain formulas for the correlation 
functions we use (|2.2|) together with our expressions for the T n . Of course the same formulas 
must result if we take the appropriate coefficients in ( p.4|) . 

It turns out that the 2x2 block matrices (^K^{y r ^ y s )^j are always self-dual. This means 
that the matrices on the diagonal are transposes of each other and the matrices off the 
diagonal are antisymmetric. Because of this the expression for the correlation function 
given by ( [2.2|) shows that (|1.2p holds when the determinant is interpreted as a quaternion 
determinant. For a discussion of this point see ||, sec. 6.2. 

4. Formulas for the spacing probabilities 

We obtain the formulas first assuming that (3 = 2 and (|3.1|) holds. To evaluate i?(0; J) 
we must integrate P^{x\, • • • , x^) over the region where all the Xi lie in the complement of 
J. Thus we simply take / = — Xj and obtain 

E(0; J) = det {I - K N Xj). 

More generally, to evaluate E(n\, ■ ■ ■ , n m ; J 1; • • ■ , J m ) we integrate P^(xi, ■ ■ ■ , x^) over the 
region where rij of the Xj lie in Jj and all the other Xj lie outside the union of the Jj. We 
see that E(ri\, ■ ■ ■ , n m ; J 1; ■ ■ • , J m ) equals the coefficient of (Ai + l)™ 1 • ■ ■ (A m + l)" m in the 
expansion about Ai = • • • = X m = — 1 of 

r- r- m 

/•••/ P N {xi, - ■ ■ ,xn) X\[l + ^2\iXjXxj)} dxx- ■ -dx N . 

j i i 
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Therefore we take f(x) = YhLx AjXj^x) in ( |3.1| ) and deduce that 

E(ni, ■ ■ ■ , n m ; Ji, • • • , J m ) 

= — ; tttt^ ^ — det (l + K N V AiXj.) • (4.1) 

n x ! • • - n m ! dA™ 1 • ■ ■ V w ^ -V A 1= ...=A m =-i v ; 

To see the connection with spacings between eigenvalues, consider the quantity 

[£(0; {x + Ax,y))-E{0; {x, y))\ - [£7(0; (x + Ax, y + Ay)) - £7(0; (x, y + Ay))], 

where Ax and Ay are positive. The difference on the left equals the probability that there 
is no eigenvalue in (x + Ax,y) but there is an eigenvalue in the larger interval (x,y), which 
is the same as saying that there is no eigenvalue in (x + Ax, y) but there is an eigenvalue 
in (x, x + Ax). Similarly the difference on the right equals the probability that there is no 
eigenvalue in (x + Ax,y + Ay) but there is an eigenvalue in (x,x + Ax). Therefore the 
difference equals the probability that there is no eigenvalue in (x + Ax, y) but there is an 
eigenvalue in both (x, x + Ax) and (y, y + Ax). Hence if we divide this difference by Ax Ay 
and take the limit as Ax, Ay — > we obtain the joint probability density for consecutive 
eigenvalues to lie in infinitesimal intervals about x and y. On the other hand this same limit 
equals 

-£r, m <*■»»■ 

Therefore this is the formula for the joint probability density, and explains why the spacing 
distributions are intimately connected with Fredholm determinants. 

To obtain the conditional probability density that, given an eigenvalue at x, the next 
eigenvalue lies in an infinitesimal neighborhood of y, we must divide the above by the prob- 
ability density that there is an eigenvalue in an infinitesimal neighborhood of x. This is 
the 1-point correlation function R\(x, x) which, we know, equals K^(x, x). Therefore this 
conditional probability density equals 

1 d 2 

-£7(0; (x,y)). 



Kn(x, x) dydx 

We now compute Pj{ ), the joint probability density that the eigenvalues con- 

tained in an interval J lie in infinitesimal intervals about x\, • • • ,x n . We observe first that 
the probability that a small interval contains more than one eigenvalue goes to zero as the 
square of the length of the interval. So we shall apply ( [4.1|) to small intervals about the 
Xi (taking these rii = 1) and the intervals constituting the rest of J (taking these = 0), 
divide by the product of the lengths of the small intervals and let their lengths tend to 0. 
We denote the small intervals by Ji,---,J n and the rest of J by J'. Denote by Ej lt ... t j n the 
probability that each Jj contains one eigenvalue and J' none. By (|4~T|) we have 

Qn n 

E Ji,-,Jn = ttt wr~ det (/ - K N Xj, +K N J2^i 

OA\ ■ ■ ■ OA n v i=1 



s 



Qn n 

We factor the operator in this expression as 

n 

{I-K N Xj) {I + Rj 5>XjJ, 

where i?j is the resolvent kernel of K~n Xj, the kernel of (/ — Xj)~ x X j. Determinants 
multiply, and the determinant of the first factor equals E(0; J), as we know. Thus 

Now we can proceed as for the correlation and cluster functions. If 

n 

A(z u ■■■,z n ) = det (1 + RjY, z i X^) ~ h 

n 

B(zi, ■■■,z n ) = log det (/ + Rj ^ z { Xj t ) , 

i=i 

then Ej lt ... j j n /E(0; J) equals the coefficient of z\---z n in the expansion of A about z\ = 
■ • ■ — z n — and A = e B — 1. The coefficients of A are determined from the coefficients of 
B by formula ( |2.4j ), which gives in this case 

As = z2B Sl - ■ ■ B Sm , 

where the sum is taken over all partitions {Si, ■ ■ ■ , S m } of S into nonempty subsets. To 
evaluate 

^ = Il7r- log det (l + RjJ2 z i x Ji) 

we use the general facts for an operator function K(z) that 



zi=—=z n =Q 



-flog det K(z) =trK(z)- 1 K'(z), ^-K(z)- 1 = -K(z)- 1 K'iz) K(z)- 1 . 
dz dz 

We apply this consecutively to the derivatives with respect to the Zi and find that if S = 
{ii, ■ ■■ ,i k } then 

B s = - — t Y tr Rj Xj. ■■■ Rj Xj- , 

where a runs over all permutations of {1, • • ■ , k}. 

Eventually we have to divide the various products Bg 1 ■ ■ ■ B$ m by the product of the 
lengths of all the Jj and take the limit as these lengths tend to 0. We may consider separately 
each Bs and its corresponding intervals since the Si are disjoint. So we compute 

tr RjXj il -"RjXj ik = / ••• / Rj{x ik ,x il )Rj{x il ,x i2 )---Rj(x ik _ 1 ,x ik )dx il ---dx ik . 
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If Ji x , • • • , Ji k are small intervals about , ■ ■ ■ , Xi k then dividing by the product of their 
lengths and letting the lenths tend to gives in the limit 

Pj (-^ife j ■Eix ) Pj ) «2 ) ■^•Z \ X ik—l ' ) " 

Now we would have obtained exactly the same results if, instead of starting with the 
determinant of the operator I + Rj J2ziXj t , differentiating with respect to the Zi, setting 
all Zi — 0, dividing and taking a limit, we had simply started with the determinant of the 
matrix 

[5ij + Rj(xi,Xj) zj) . . 

and differentiated with respect to the Zj and set all Zi = 0. It follows that the the ratio 
Pj(x u ■ ■ -,x n )/E(0; J) equals 

— det (Sij + Rj(xi, Xj) zj) 

dz\---dz n v J *i=-=* n =o 

Therefore we have obtained the formula 

Pj(xi, ■■■,x n ) — E(0; J) det (Rjfa, Xj)) . . . (4.2) 

Note the similarity to formula (|1.2| ) for the correlation functions: to obtain the formula for 
the ratio Pj(xi, ■ • ■ , x n )/i?(0; J) we replace the kernel K N (x, y) by the kernel Rj(x, y). 

We mention that rather more complicated formulas for these probability densities were 
obtained by Mehta and des Cloizeaux [f|. (See also ||.) These involve the eigenvalues and 
eigenfunctions of the operator Xj but they must must be equivalent to ( |4.2j ). 

For the spacing probabilities when (3 = 1 or 4 we must make modifications as at the end 
of the previous section. We have in this case 

E(0; J) 2 = det (I - K N Xj) 

and the ratio Pj(xi, ■ ■ ■ , x n )/E(0; J) is again given by the same formulas as for the cor- 
relation functions but with K^(x,y) replaced by Rj(x,y). One can show that the block 
matrices (Rj(xi,Xj)) that arise here are also self-dual, and this implies that the formula for 
Pj(xi, ■ • • , x n ) is as above if the determinant is interpreted as a quaternion determinant. 



5. The (3 = 2 matrix ensembles 



Consider the Hermitian case first. We use ( |1.3| ) to express 

/" ' / ~ x kf W. w { x i) Wi l + f( x j)) dx i ■ "dx N (5.1) 



determinant. Since 



j<k 



Y\_{ x j — x k) — det (x{)j = o t ... t N_i t k=l,-,N 
3<k 
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applying (|1.3|) with (pj(x) = i/ij(x) = x J and with the measure w(x) (1 + f(x)) dx shows that 
this equals N\ times 



det ( f x j+k w(x) (l + f(x))dx 



j,k=0,-,N-l 



If we denote by {(fj(x)} the sequence obtined by orthonormalizing the sequence {x 3 w(x)^} 
then we see that the above is equal, except for a different constant factor depending only on 
N, to 

det ( / (fj{x) tp k (x) (l + f(x)) dx) 



>j,k=0,-,N-l 

= det (d jik + J ipj(x) if k (x) f(x) dx ) jk=0 .,. jN _ 1 - ( 5 - 2 ) 

Now we use something which is needlessly fancy but which is very useful, namely the 
general relation det (/ + AB) = det (/ + BA) for arbitrary Hilbert-Schmidt operators A and 
B. They may act between different spaces as long as the products make sense. In our case 
we take A to be the operator from L 2 to with kernel A(j,x) = <fj(x) f(x) and B the 
operator from to L 2 with kernel B(x,j) = <fj(x). Then 



AB(j, k)= tpj(x) ip k (x) f(x) dx, 
so (|5.2j ) equals det (J + AB). Now BA is the operator on L 2 with kernel K^(x, y) f(y) where 

N-l 

K N (x,y) = Vk(x)ip k {y). 

Therefore ( ^.2|) also equals det(J + K N f). 

Hence we have shown for these ensembles that (|3.1|) holds up to a constant factor de- 
pending only on N. In fact the constant factor now must be 1, as is seen by setting / = 0, 
and so ( |3.1|) is true as it stands. 

It is easily seen directly that the nonzero eigenvalues of K^f are the same as those of 
the matrix (/ (pj(x) <p k (x) f(x) dx), for any eigenfunction must be a linear combination of the 
ipj and solving for the coefficients leads directly to the matrix. So the fact det(J + AB) = 
det(J + BA) is very simple here. 

When the measure f(x) dx is discrete, in other words when f(x) = J2r=i z r $( x ~Vr), the 
integral in the right side of ( |5.2| ) is replaced by J2r=i z r fjiVr) Vk(y r )- If we set 

A(j,r) = z r (pj(y r ), B{r,j) = (pj{y r ), (j = 0, • • • , JV - 1, r = l,---,n) 

then our matrix is AB, and 

BA(r, s) = K N (y r , y s ) z s , (r, s = 1, - • • , n) 



with K N as before. Thus in ( |3.1| ) the operator / + K N f is replaced by the matrix (|3.2|), as 
claimed. (The same argument will hold for the other ensembles, and there will be no need 
to repeat it.) 



11 



For the circular (3 = 2 ensemble ( |5.1| ) is replaced by 

/•••/ 11^ '''' ~< '''') T[(e iX3 ~e iXk ) '[[{l + f{x j ))dx 1 ---dx N . 

j<k j<k j 

In the application of (|1.3| ) we take 4>j (x) = e~ lXj / V2~rr, ipj (x) = e %Xi j \f2~n so that the analogue 
of the matrix in (15. 21) becomes 



— / e l{k ~ j)x f(x)dx. 
(The 1/27T factor is needed to obtain when / = 0.) Now we set 

A(j,x) = l-e-^-^fix), B(x,j) = e^-W-, 

(the extra exponential factors yield a simpler formula), so that the above matrix equals AB 
while BA equals the operator with kernel K^{x,y) f(y) where now 

K N (x y) = IV>¥)M = J- Sin f (x ~^ 
K N [x,y) 2n ^e 27rsin|(x-y)- 

Thus (|3.1|) holds for the circular ensemble with this replacement for K^. 



6. The (3 = 4 circular ensemble 

We begin with the evaluation of 

j ... [ TJ|e^-e^| 4 n(l + /(^))^i---^iv (6.1) 

Pfaffian. It is a pretty fact that 

T[(xj - x k f = det (4 jxl 1 ) . (6.2) 

11 v J ' \ k J k ) j=Q 2N-1, k=l,-,N v ' 

j<k 

This is seen by writing the product representation of the Vandermonde determinant 

det (xi yi] 

V K aK ) j=0,-,2JV-l, k=l,-,N 

then differentiating with respect to each y k and setting y k = x k . 
If we replace each x k by e tXk and use the fact that 

| e i% _ e te fc | = ^ e -i=4^ ( e »i _ e iXk ) (6.3) 
if < x k (both lying in the same interval of length 27r) we see that 

TJ \e ix i — e ix fc| 4 = e -2i(JV-i) xj TT^g^j _ e ix ^ . 

j<k j<k 
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By ( |6.2|) this is equal to 

e -MN-DE*i det (e ijx * je^ = det (e^-^+l)** 

V J I j=0,...,2JV-l, k=l,-,N y > 

= det (e^'-^+s)** (j-N+ i) e ^- N+ ^ x A = det (e* p:Cfc pe ijxCfc ), 

where in the last determinant = 1, • • • , N as before but p runs through the half integers 
-N+l-N+l---,N-±. 

By the symmetry of the integrand in (|6.1| ) the integral equals TV! times the integral over 
this region x\ < • • • < xn ■ So we can use formula ([Ol ) to deduce that the square of (|Q| ) 
equals a constant depending only on N times the 2N x 2N determinant 

det (f(q-p) e l{p+q)x (l+f(x))dx). 



Both indices p and q run over the half- integers —N + |, • • • , N — |. If we reverse the order 
of the rows and divide each column by its index q we see that this determinant is equal to 
another constant depending only on N times 

d6t {^J( 1 + l) el{ ~ P+9)X ( 1 + f( x )) dx )= det (^ + ^ji l + ^)^' P+q)X f^) d ^)- 

Here we see the sum of products e~ tpx e iqx + p e~ ipx - e tqx referred to in the introduction, 
and we write it as a matrix product 



( e~ ipx ipe~ ipx ) y f Q J qx j . 
(The reason for the insertions of i will become apparent.) Thus if we set 

e iqx 

iq 



A(p,x) = —f(x) ( e %px ipe vpx ) , B(x,q) = ( f 



then the above matrix is I + AB. In this case BA is the integral operator with matrix kernel 
K N (x,y) f(y) where 



K N (x,y) 

If we write 



/ —J2 e*^ - ^ e ip( - x ~ y ^ 

J_ J_ p ip(x-y) J_ ip(x-y) 
\ 4tt ^ ip ° 471^° 



1 1 sin Nt, d r x 

Sn{x) = — £ e ipx = —— -T-, DS N (x) = —S N (x), IS N (x) = / S N (y) dy, 

Z7T p 2lX Sill ±X dx JO 

then 



-,/ S N (x-y) DS N (x-y) 
K N (x,y) = -\ | • (6.1) 

\ IS N (x - y) S N (x - y) 
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Note the relationship between Sn(x — y) and the kernel K^{x ) y) which arises in the (3 = 2 

ensemble: the former equals the latter with N replaced by 2N . 

Recall that we have been computing, not the integral ( |6.1|) , but its square. Thus we have 
shown for this ensemble and with this that ( |3.3| ) holds up to a constant factor depending 
only on N. As before, taking / = shows that the constant factor equals 1. 



7. The (3 = 1 circular ensemble 



We assume here that N is even and begin with the evaluation of 

/■■•/ II ^< l[{l + f(x j ))dx 1 ---dx N . (7.1; 

j<k j 

Using (|6.3| ) we see that in the region x± < ■ ■ • < xjv, 

TJ \e ix i — e iXk \ = j-NiN- 1 ) e ~ iIi f 1J l2 x j Y\_{e iXk — e iXj ) 



-N(N- 



- 1 ) e -i V E x i det f e ijXk ) = i-W- 1 ) det f e ipXk ) . 

V / jf=0,---,JV— 1, k=l,-,N V /' 



N i 1 



where in the last determinant k = 1, ■ ■ ■ , N and p runs through the half integers from — y + 
to y — |, as before but with N replaced by N/2. The z-factor in the last term equals ±1. 



We now apply ( |1.4|) with the measure (1 + f(x)) dx, using the fact that n \ e%Xi ~ e %Xk \ is a 
symmetric function, and conclude that the square of ( [7.1|) is a constant depending only on 
N times 

det( / f sgn(y-x)e ipx e^(l + f(x))(l + f(y))dydx). 



First, we replace p by — p as before. Then we set e(x) = \ sgn(x), and replace sgn(y — x) by 
e(x — y) in the determinant. These things just change the multiplying constant. We shall also 
use the notation e p (x) for the normalized exponential e ipx /y/2n. With these modifications 
the above becomes 

det( J J e(x-y) e_ p (x) e q (y) (1 + f{x)) (1 + f(y)) dydx). 

Denote by e the operator with kernel e(x — y). Observe that e is antisymmetric and 
ee p = e p /ip. (The latter uses the fact that p is an integer plus |.) Writing 

(1 + f(x)) (1 + f(y)) = 1 + f(x) + f(y) + f(x) f(y) 

and using the fact that e is antisymmetric, we see that the double integral equals 

1 r 1 1 

—8 P , q + / [— f(x) e_ p (x) e q (x) dx + — f(x) e_ p (x) e q (x) - f(x) e q (x) e(fe^ p )(x)\ dx, 
iq J iq ip 

so the determinant equals another constant depending on N times that of the matrix with 
p, q entry 

$P,q + / lf( x ) e -p( x ) e q( x ) + ~ f( x ) e- p (x) e q (x) - iq f(x) e q (x) e(fe- p ){x)\ dx. 
J p 
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At this point setting / = shows that the multiplying constant equals 1. 
We write the integrand in the p, q entry of the matrix as a matrix product 



( k f( x ^ e ~P^ ~ f( x ^ £ (f e -p)( x ) f( x ) e -p( x ) ) 



iq e q [x) 



e q (x) 



and apply / + AB — ► I + BA as before to replace this by a matrix kernel. For simplicity of 
notation we shall write out the operator rather than the kernel and use the notation a <g> b 
for the operator with kernel a(x) b(y). The operator / + BA is equal to 

/ / + £ e p <g> /e_ p - E *P ep ® / e(/e_ p ) E *P e p <g> /e_ p 

V E 5 e p ® /e_ p - E e p <g> / e(/e_ p ) / + E e p ® /e_ p . 

Now it is an easy fact that for an operator A we have (a ® 6) A = a ® (A* 6). Therefore 
the second sums in the entries of the first column of the matrix, with their minus signs, equal 
the corresponding entries of the second column right-multiplied by the operator ef. (Here 
/ denotes multiplication by the function, which is symmetric, and we used the fact that e is 
antisymmetric.) It follows that the matrix may be written as the product 

/ I + E e p <S> /e_ p E ip e P ® /e_ p 

\Ei e p® f e -p ~ £ f I + J2 e p ® fe- p 

The operator on the right has determinant 1 and so the determinant of the product equals 
the determinant of the operator on the left, i.e., the determinant of 




p Ky d—p 



E ipe 



p vv C— p 



E I 



6_ p £ E ^ 



p C_ p 



If we recall the definition of e p and the notation of the last section we see that (|3.3| ) holds 
for this ensemble when has kernel 



K N {x,y) 



SN_(x-y) 



DSn(x —y)\ 



v ISn(x - y) - e(x - y) Sir(x-y). ) 



8. The (3 = 4 Hermitian ensembles 



We begin now with the iV-fold integral 

/"'/ Ii( x 3 ~ x kY II u; ( :r i) IK 1 + f( x j)) dx i ' dx 



j<k 
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Using Q5.2| ) as it stands and ( |1.5| ) we see that the square of this iV-fold integral equals a 
constant depending only on N times 



det (J{k- j) x j+k ~ l (1 + f(x)) w(x) dx) 



j,k=0,---,2N-l 



We replace the sequence by any sequence {pj{x)} of polynomials of exact degree j. 
Except for another constant factor depending only on N, the above equals 

det ( J (pj(x)p k (x) - p'j(x) p k (x)) (1 + f(x))w{x) dx). 

Just as when /3 = 2 it was convenient to introduce functions ipj which were equal to the 
orthonormal polynomials times so now we introduce ipj = Pjiu? . The matrix is equal to 



(ipj (x)ip' k {x)- ipj (x) ip k (x))(l + f (x) ) dx) , 
the extra terms arising from the differentiations having cancelled. We write this as 



M + [ J (iPj(x) iP' k (x) - iP'j(x) ip k {x)) f(x) dx 

where M is the matrix of integrals f(ipj(x) ip' k (x) — ip'j(x) ipk(x)) dx. 

Next we factor out M, say on the left. Its determinant is just another constant depending 
only on N. If M _1 = (fJ.jk) and we set f]j{x) = E Hjkipkix) then the resulting matrix is 



1 + ( / (Vj( x ) ^k( x ) ~ Vj( x ) i/>k(x)) f( x ) dx), 
Next, we write the sum of products f]j(x) ip' k (x) — t]j(x) ip k ( X^j clS db matrix product 

V ^Pk(x) 



T]j(x) -ij{x) 



Thus if we set 



= /0) ( r]j(x) -ij{x) ) , B(x, j) 



/ iP'j(x) 



then the above matrix is / + AB. In this case BA is the integral operator with matrix kernel 
K N {x,y) f(y) where 



K N {x,y) 



I E^'j(x)vj(y) -E*P'j(x)rij(y) 



\ E *Pj 0) Vj (y) - E ipj ( x ) v'j (y) 

The sums here are over j = 0, • • ■ , 2N — 1. 
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We have shown that with this Kn ( |3.3| ) holds up to a constant factor depending on N 
and as usual taking / = shows that the constant factor equals 1. 

Recall the definition of the rjj. If we write 

2N-1 

S N (x,y)=2 J2 *P'j( x ) Vjk*Pk(y) 

j,k=0 



and 



then 



2N-1 2N-1 

IS N (x,y) = 2 ^A x ) Mki>k{y), S N D(x,y) = -2 ^j( x ) VjktfM 
j,k=o j,k=o 



K N {x,y) = - 



( S N {x,y) S N D(x,y) 



i.l) 



\IS N (x,y) S N (y,x) 

The explanation for the notation is that if Sn is the operator with kernel SN{x,y) then 
SnD(x,v) is the kernel of SnD {D = differentiation) and ISn(x,v) is the kernel of ISn 
(I = integration, more or less). The factors 2 and \ are inserted to maintain the analogy 
with 0. 

One might wonder about the choice of the pj and the resulting fijk- Since M is anti- 
symmetric the formulas would be simplest if it were the direct sum of N copies of the 2x2 

matrix Z = j),who S e inverse equals its negative. The polynomials that achieve 

this are called skew-orthogonal. But there is no actual necessity for using these since, as is 
easily seen, any family of polynomials leads to the same matrix kernel. 



9. The (3 = 1 Hermitian ensembles 



Again we assume that iV is even and begin now with 

/■ "/ n \ x j - x k\ n^fe) IK 1 + f( x j)) dx ^ ~" dx 



This is equal to iV! times the integral over the region x\ < ■ ■ ■ < xn and each \xj — Xk\ = 
Xk — Xj there. Using ( |1.4|) we see therefore that the square of the integral equals a constant 
depending only on iV times 



det ( / / e(x - y) x 3 y k (1 + f{x)) (1 + f{y)) w{x) w(y) dydx 



j,k=0,--,N-l 



As in the last section we replace by an arbitrary sequence of polynomials {pj(x)} of 
exact degree j.. The above becomes a constant depending only on times 

det( / / e(x-y)pj(x)p k (y) (l + f(x) + f(y) + f(x) f(y)) w(x)w(y) dydx). 
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Now set fa(x) = Pj(x) w(x) and denote by M the matrix with j, k entry 

J J e(x-y)p j (x) p k (y)w(x)w(y)dydx. 
Then the last determinant equals 

det (M + y Je(x-y)fa(x)fa(y)(l + f(x) + f(y) + f(x)f(y))dydx) 

= det (M + J [f fa ei) k - ffa efa -ffa e(ffa)} dx) . 

All expressions in the integrand are functions of x. As earlier, the operator with kernel 
e(x — y) is denoted by e. 

As in the last section, we factor out M on the left, write M _1 = (//_,-&) and set rjj = 
E Hjkipk- What results is the determinant of 



/ + 



( J [f Vj eipk - f4>k erjj -ffa e(fr)j)] dx) , 



and, as usual, at this point the multiplying constant equals 1. 

The integral is equal to the j, k entry of the product AB, where 

A(i,x) = ( -ferjj - fe(frij) fvj ) , B(x,j) = 



fa 
efa 



The operator / + BA equals 

/ I -Efa ® ferij -Efa ® / s(fVj) Efa®fVj \ 

V - E efa ® / eVj ~ E efa ® / e(/?fr) I + Eefa®f Vj J 
As in the circular ensemble, we can write this as the matrix product 

I-Efa^ferij T,fa® f Vj \ ( 1 

-J2efa ®fer }j - ef I + J2efa ®fVj) \ ef 1 

The determinant of the product equals the determinant of the first factor, i.e., the determi- 
nant of 

/ - E fa <8> £Vj E fa <8> Vj \ 

i+ /■ 

V - E efa ® £Vj ~ e E efa <8> Vj ) 
Recalling the definition of the qj we now define 

N-l 

S N (x,y) = - J2 fa( x ) Pjk&PM, 

j,k=0 
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JV— 1 N-1 

IS N (x,y) = - ^j(,x)njkeil>k(y), S N D(x,y) = *Pj( x ) Vjk ^k{y), 

j,k=0 j,k=0 

and we see that the (|3.3|) holds for this ensemble when Kn has kernel 



K N {x,y) 



S N (x,y) S N D(x,y) 
\ IS N (x,y) - e(x - y) S N (y,x) 



(9.1) 



Once again the matrix kernel is independent of the choice of the pj but the expression 
is simplest when the polynomials are skew-orthogonal in this context, i.e., when M is the 
direct sum of N/2 copies of Z. 



10. The Gaussian ensembles 



The Gaussian symplectic ensemble ((3 = A). The weight function w(x) now equals e~ x . 
As in sec. 5 we use the notation ipj for the polynomials orthonormal with respect to the 
weight function times the square root of the weight function. It turns out that in this case 
the skew-orthogonal polynomials are simply expressed in terms of the <pj. In fact, in the 
notation of sec. 8, we have 

^2n+l = f2n+l/ V2, 1p 2 n = ~^2n+l/v / 2- 

In other words, 

P2n+l(x) = e x2/2 ip 2 n+l(x)/V2, P 2 n(x) = -e^ 12 £<p 2n +l/ V5. 

The point is that p2 n as so defined is actually a polynomial of degree 2n, as is easily seen. 
The skew-orthogonality is easy. If j and k are of the same parity the corresponding matrix 
entry equals 0, and if they have opposite parity we compute, (using the fact that De = 
Identity) 

J (P2n(x) p' 2m +l( X ) ~ P2n( X ) P2m+l(x)) dx = J (tfj 2n 1p 2m +l ~ V4n ^2m+l) dx 

= \ J (</?2n+l V2m+1 ~ £^2n+l <f' 2m +l) d X. 

Integration by parts applied to the second integrand and the orthonormality of the (pj show 
that the above equals b n . m . 

We now compute the entries of the matrix ( |3 . 1|) . If we keep in mind that the inverse of 
M in this case equals — M we see that Sn(x,v) is equal to 1/2 times 

N-1 N-1 

<f2n+l{x) (P2n+l(y) ~ ^ ^2n+l{ X ) ^2n+l(y)- (10.1) 
n=0 n=0 

To compute the second sum we use the differentiation formula 
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If (ciij) is the antisymmetric tridiagonal matrix with a^j-i = yi/2 then 

f'i = °ij <Pj> <Pj = a Ji £l Pi- ( 10 - 2 ) 
j>0 i>0 

The first is just a restatement of the differentiation formula and the second follows from the 
first by applying e and interchanging i and j. By the first part of (|10.2| ) the second sum 
in ( |10.1| ) equals J2 a ij <Pj ' £ Vi summed over all i, j > such that i is odd and < 2N — 1. 
Using the antisymmetry of (a^) and a different description of the range of summation, we 
see that this equals — Yjfj ' a ji e{ Pi summed over alH, j > such that j is even and < 2iV 
except for the single term corresponding to i — 2N + 1, j = 2N. (Recall that = unless 
\i — j\ = 1.) By the second part of ( |10.2| ) this equals 

N 



n=0 



Recall that Sjy(x,y) equals 1/2 times the sum (|10.1| ). If we now denote by S^f(x,y) the 
sum itself, then 



2N 

S N (x,y) = J2 ¥>n{x)Vn{y) + \JN +i (pwix) eip 2N+1 (y), 

n=0 

and consequently 



2iV 



IS N (x,y) = £<fn(x) ip n (y) + y N +± £^ 2 7v(^) £<£> 2 iv+i(2/), 

n=0 



2iV 



S N D(x,y) = -J2 <Pn(?) <p' n (y) - -y/W+i ^jvfa) ^jyv+ifa). 

n=0 

(That ISN(x,y) is obtained by applying e to ( |10.1|) as a function of x is a consequence of 
the fact that only odd indices occur there.) Hence our matrix kernel is given by 



1 f S N (x,y) S N D(x,y) 

K N {x,y) 

IS N (x,y) S N (y,x) 
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in analogy with 



The Gaussian orthogonal ensemble ((3 = 1). In order to continue using the same functions 
(fj as before we shall now take as our weight function w(x) = e~ x I 2 . Thus the functions ipj 
of sec. 9 will be of the same form, polynomial times e~ x I 2 . For the pj to be skew-orthogonal 
in the present context the antisymmetric matrix M with j, k entry 

e(x - y) ipj(x) ipk(y) dydx = J ipj{x) eip k (x) dx 

must be a direct sum of N/2 copies of Z. (Recall that N is even.) The ipj are again very 
simple. They are given by 

1p2n = V?2n, 1p2n+l = ^' 2n - 
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That these give polynomials of the right degree (and parity) is clear, and we compute 
The sum S^(x,y) in ( |9.1|) is now given by 



2 1 



Y (lp2n(x) £lf>2n+l(y) ~ ^2n+l(x) #2n(?/)) = Y ^2n{x) <P2n(v) ~ Y ^2n( X ) ^2n{y) ■ 
n=0 n=0 n=0 



By ( |10.2j ) the second sum on the right, equals a ij ¥>j • £( Pi summed over all i, j > such 



that i is even and < N — 2. This equals — fj • a ji £L P% summed over all i,j > such that 
j is odd and < N — 1 except for the single term corresponding to i = N, j = N — 1. This 
sum equals 



JV 
2 



N 

Y <P2n-l(x) <p 2 n-l(y) ~ \ Ht <PN-l( X ) £l PN{y)- 
n=0 V A 



Hence we have now 



N-l 



SN(x,y) = Y Vn{x)Vn{y) + \ -~ <fiN-l(x) S(f N (y), 
and consequently the other entries of the matrix K^{x,y) in (|9.1|) are given by 



N-l 



IS N (x,y) = Y £<Pn(x) <fn{y) + J-?r £(fN-i{x) e<p N (y), 

n=0 » 1 



N-l 



1 N 



S N D(x,y) = - Y <Pnfa) <p' n (y) - \ — <p N -i(x) ip N (y). 

n=0 » / 
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